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INTRODUCTION 


The increasing interest in minimum weight designs lor aeronautical and aerospace structures has 
generated substantial interest in the analysis of the elastic stability and postbuckling behavior of 
structures subjected to inplane compressive loads. For thin homogeneous plates, classical plate theory 
predicts deformations and inplane stresses that are comparable to those of three-dimensional elasticity. 
Transverse stresses in thin plates are generally small compared to inplane stresses, and thus, both 
classical theory and first-order shear deformation theory give satisfactory results. However, since both 
theories are two-dimensional, they are not accurate enough to predict transverse stresses directly. 
Accurate nonlinear theories are required for the analysis of thick plates in which these transverse 
stresses become more significant. 

It is often sufficient to use an accurate nonlinear two-dimensional theory to solve some three-dimensional 
nonlinear elasticity problems. However, when through-the-thickness effects become more dominant, it 
is important to use a nonlinear theory that takes into account such effects. One such theory has been 
derived in reference 1 for laminated and thick plates with three-dimensional flexibility effects. This theory 
can predict directly the transverse stresses as well as the inplane stresses by using trigonometric terms in 
addition to the usual constant and linear terms representing through-the-thickness variation of the 
displacements. However, this theory cannot satisfy the surface boundary conditions of a plate. 

The purpose of the present paper is to present the results of an investigation of the buckling and 
postbuckling response of orthotropic plates loaded in compression. Classical nonlinear von Karman 
theory using the Kirchhoff assumptions and three nonlinear transverse shearing theories are used to 
predict results for different values of plate width-to-thickness ratios in the postbuckling range. The 
nonlinear transverse shearing theories are: first-order shear deformation theory, references 2 and 3, 
higher-order shear deformation theory, reference 4; and three-dimensional flexibility theory, reference 
1 . The idea of satisfying exactly the static tangential or surface boundary conditions on the external 
planes of the plate (or shell) was used in references 5, 6, and 7. The first papers dealing with 
postbuckling where the static tangential or surface boundary conditions on the external planes of the 
plate were satisfied are references 8 and 9. The present derivation of the higher-order shear 
deformation theory has the advantage of having nonlinear through-the-thickness terms without 
contributing additional unknowns to the first-order shear deformation theory. In addition, it satisfies the 
surface boundary conditions of the plate. The essential difference between the higher-order shear 
deformation theory and the three-dimensional flexibility theory is that the higher-order shear deformation 
theory is a two-dimensional theory that uses cubic terms, whereas the three-dimensional flexibility theory 
is a three-dimensional theory that uses trigonometric terms in addition to the constant and linear terms 
that represent the through-the-thickness variation of the inplane displacements. The paper presents the 
derivation of the nonlinear plate equations for buckling of plates loaded in axial compression for both 
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higher-order theories. This paper also presents postbuckling results for the average longitudinal 
compressive direct stress resultant and maximum stress resultants as a function of the applied 
displacements, and maximum out-of-plane displacement as a function of the applied end-shortening 
displacement. The plates considered in this paper are infinitely long with side edges simply supported 
and are loaded in uniaxial compressive end shortening. The side edges are free to slide along the edges 
to allow constant longitudinal strain. Results for the four theories are presented for aluminum plates and 
for composite plates with a symmetric lay-up composed of many layers of thin ±45° plies. 
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SYMBOLS 

Plate extensional stiffnesses 

Dimensions of the rectangular plate parallel 
to x, y, and z axes, respectively 

Stiffnesses used in Hooke's Law 

Plate bending stiffnesses 
Young's modulus 

Plate stiffness components, defined in equation 
(A14) in Appendix 

Moment resultants in the plate 

Inplane stress resultants in the plate 

Transverse stress resultants in the plate 

Average compressive load per unit length 

Value of N xav at buckling 

Functions of y defined in the Appendix 

Applied end shortening 

Value of U at buckling 

Displacements in x, y, and z directions, 
respectively 

Plate coordinates 

Functions of y defined in the Appendix 
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Strains in the plate 
Buckle half-wavelength 
Poisson's ratio 
Stresses in the plate 

THEORY 

A brief outline of the derivation of the four different theories compared in this paper is presented in 
this section. The derivation of equations using classical von Karman-Kirchhoff theory has been 
presented in reference 9. The derivation of equations using first-order shear deformation theory has 
been presented in references 2 and 3. The derivation of the equations for the two higher-order 
theories are not given in detail elsewhere, so they are presented in the appendix. The general 
approach used in deriving the equations to be solved is the same as in reference 10. First, the 
displacement functions for each theory are identified. Then the nonlinear strain-displacement 
relations are written to include the assumption that the displacements are sinusoidal along the length 
of the infinitely long plate. Stress-strain relations are defined for a "specially orthotropic" plate. 
Application of the principle of virtual work leads to ordinary differential equations and variationally 
consistent boundary conditions which are solved using a procedure based on Newton's method as 
discussed in reference 1 1 . 

The displacements considered for each theory are: 

Classical von Karman-Ki rchhoff theory 

u(x,y,z) = u°(x,y)-w,J(xy)i 
v(x,y,z) = v°(x,y) - w,y(x,y) jjj- 
w(x,y,z) = w°(x,y) 

First-order shear deformation theory 

u(x,y,z) = u°(x,y) + u Vy) 
v(x,y,z) = v°(x,y) + vVy)^ 

w(x,y,z) = w°(x,y) 


e x . e y- e z> Yyz> Txz' Txy 

X 

P 

CT X , <Ty. &Z’ T yz- X XZ’ T xy 


( 2 ) 



Hiaher-order shear deformation theory 


u(x,y,z) = u (x,y) + 


u Vy)-y- +w °(x,y)j (i) 


z_ 

h 


v(x,y,z) = v (x,y) + 


w(x,y,z) = w°(x,y) 


'Vy)-^ 5 ^^ + w .y( x .y)j j^-J 


h 


( 3 ) 


Three-dimensional flexibility theory 

u(x,y,z) = u°(x,y) + u a (x,y) ^ + u s (x,y) sin £2. 

h h 

v(x,y,z) = v°(x,y) + v a (x,y) -2- + v s (x,y) sin 22- (4) 

h h 

w(x,y,z) = w°(x,y) + w c (x,y) cos^ 

h 

In this paper the zero superscripts correspond to the constant-in-z terms, the a superscripts 
correspond to the algebraic-in-z terms, and the s and c superscripts correspond to the 
trigonometric-in-z terms. 

Both the classical von Karman-Kirchhoff and the first-order shear deformation theories have inplane 
deformations u and v which are linear in z. Classical theory, however, has the additional assumption 
that there is zero transverse shearing (y xz = = 0), thus eliminating u a and v a in favor of 

derivatives of w°. 

The higher-order shear deformation theory considers inplane deformations u and v which are cubic 
in z. As explained in reference 10, the squared-in-z term vanishes and the cubic term does not 
introduce any new variables beyond those that appear in first-order shear deformation theory if the 
boundary conditions are satisfied at z = ± h/2. The three-dimensional flexibility theory considers 
trigonometric terms in u, v, and w beyond the expressions considered for the deformations of first- 
order shear deformation theory. 

To account for the applied displacement U, 

u°(x,y) = - U - + U2 (x,y) 
a 
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v°(x,y) = v 0 (y) +V 2 (x,y) 


( 5 ) 


where numbered subscipls for loads, displacements, and curvatures indicate a y dependence only. 
To satisfy the assumption that the displacements are sinusoidal along the length 

0 0 O iry 

U 2 = u 2 (y)sin^- 

X ( 6 ) 

0 0 2 jcx 

V 2 = v 2 (y)cos-^— 


All the other u coefficients can be expressed as functions of y multiplied by cos nxfk, where X is 
the half-wavelength of the buckled plate. All the other v and w coefficients can be expressed as 
functions of y multiplied by sin ioUX. The strain-displacement relations used are 


1 2 
£x=U'X + 2 W 'x 

1 2 

e y =v, y + 2’ w, y 

e z = w, z 
Yyz =v 'Z + w -y 
Txz= u -z + w -x 
Txy =u -y +v 'X + w -x w, y 


( 7 ) 


Hooke's law that relates stresses to strains for a "specially orthotropic" plate is used here 
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( 8 ) 


Ordinary differential equations and variationally consistent boundary conditions are derived using the 
principle of virtual work, and the equations are solved by Newton's method. The principle of virtual 
work applied to the internal forces of a three-dimensional body considered here is 



(cr x 8e x + Oy5ey + a z 8e z + Tyz$7yz + t x Z 8Yxz + 'c X y8y X y)dzdy dx 


( 9 ) 
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with simple support boundary conditions at y = 0 and y = b. 


The half-wavelength X of the assumed deformations for the infinitely long plates considered is 
chosen to minimize the buckling load for each given applied deformation. 

The principle of virtual work requires that the geometric boundary conditions be satisfied. Including 
additional terms in the representation of the through-the-thickness variation of the inplane 
displacements will lead to convergence and satisfaction of natural boundary conditions in the limit if a 
complete set of terms is used. An alternate approach is to use terms that satisfy the natural boundary 
conditions directly. A complete set of these terms also leads to convergence. For the present problem, 
the three-dimensional flexibility theory uses terms that do not satisfy the natural boundary conditions. 

For the higher-order shear deformation theory, coefficients of u and v in the assumed displacements of 
equation (3) are chosen such that Txz = 0 and Yyz = °- 1116 coefficients are written in terms of the 

existing unknowns u a , v a , and w° in a form which satisfies the natural boundary conditions at the top 
and bottom surfaces of the plate. Comparisons of results are valid whether or not natural boundary 
conditions are satisfied. 


RESULTS AND DISCUSSION 

The results obtained in this study for long aluminum plates with the geometry shown in figure 1 are based 
on values of Young’s modulus E = 10.7 x10 6 psi and Poisson’s ratio p = .33. The results obtained in 
this study for infinitely long composite plates with a symmetric lay-up composed of many layers of thin 
±45° plies are based on the stiffness properties 

A 11= A 22= 620340 lb/in - A 1 2= 446060 Ib/in. A^ 59000 Ib/in. 

a 44= a 55= 50000 lb/in A g6= 483520 Ib/in. 

D 11= D 22= 518 - 6 lbin Di 2 = 372.91 lb in. D00= 404.23 lb-in. 

for h = .1 in. and for any value of b. 

Buckling results given by the four theories are presented in figure 2 for finite aluminum plates and the 

results show the variation of the buckling stress coefficient with width-to-thickness ratio b/h for a range of 
length-to-width ratios a/b. Asymptotes to the curves in figure 2 give N xcr for plates of infinite length. 

The differences in the buckling results for aluminum plates with width-to-thickness ratios less than ten, 
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illustrate the need for including the effects of transverse shear deformations when determining the 
compressive buckling stress of these plates. 

Results in the form of average axial stress resultant Nx av versus applied end-shortening U, for the four 

theories and different values of the width-to-thickness ratio b/h, are presented in figures 3a and 3b for 
aluminum plates and composite plates, respectively. In figure 3a, only one curve is shown for a b/h value 
of 100 since the corresponding results for each theory are approximately the same. Even for the thicker 
aluminum plates there are only slight differences in the results given by the different theories. The 
results presented in figure 3b illustrate the nature of the more compliant ±45° composite plate. The more 
pronounced separation in the postbuckling branches of the curves for a given value of the b/h ratio 
indicates the lower transverse shear stiffness in the composite plates. 

The results show that the higher-order shear deformation theory gives a better approximation of the 
effects of shear deformation in thick plates than the first-order theory, but it is still a two-dimensional 
theory. Significant improvement beyond that already obtained with the higher-order theory requires a 
three-dimensional theory. The three-dimensional flexibility theory is an attempt to meet this requirement. 
However, the present formulation of this three-dimensional theory has limitations associated with the 
assumptions made on the w displacement of equation (A10). Although it does include a trigonometric 
term in z, this term does not satisfy the bounding conditions at the upper and lower surfaces of the 
plates, and therefore does not allow a z to have its proper influence on w. Additionally, the assumption 

that the nonlinear terms involving w c could be neglected because they were small in comparison to 

similar terms involving only w° may be questionable. A secondary effect is the choice of the specially 
orthotropic material which foregoes the influence of e z on Oy. Neither of these two effects becomes 

significant until the effects of shear deformation become more dominant. Evidence of these limitations 
are especially noticeable for results of the three-dimensional flexibility theory at a b/h value of 20 in 
figures 4b, 5b, and 6b. 

The results presented in figures 4a and 4b for the aluminum and composite plates show that the higher- 
order shear deformation theory gives the lowest value of normalized compressive Nymax for hi 9 her 
values of normalized end-shortening U . In figure 4b, the results for the three-dimensional flexibility 
theory at a b/h value of 20 exhibit nonlinear behavior in the postbuckling range. These results suggest 
an increase in the importance of the unsatisfied bounding conditions at the upper and lower surfaces of 
the plate. Similar behavior for the three-dimensional flexibility theory at a b/h value of 20 is shown in 
figures 5b and 6b. 


The results presented in figure 5a show that the normalized Nxymax as a function of normalized end- 
shortening U is nearly independent of the width-to-thickness ratio b/h for the aluminum plates, whereas 
figure 5b shows that the effect of the width-to-thickness ratio b/h is more significant for the composite 
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plates. Results for the normalized maximum deflection wmax presented as a function of normalized end- 
shortening U in figures 6a and 6b, for aluminum plates and composite plates respectively, show that the 
value of the deflection becomes increasingly dependent upon the width-to-thickness ratio b/h as the 
value of the normalized end-shortening U increases. These results indicate that shear deformation 
effects are more dominant for the composite plates, especially at lower values of the width-to-thickness 
ratio b/h. 

Present results indicate that three-dimensional flexibility theory gives lower buckling loads than the other 
theories, and produces acceptable results except when the effect of the missing nonlinear terms 
involving w c and the influence of the bounding conditions become dominant. The three-dimensional 
flexibility theory has the potential for permitting the development of a rigorous approach for obtaining 
direct through-the-thickness stress components without the current limitations caused by using 
additional trigonometric terms in z in the expansion of the transverse displacement w and by retaining 
the currently neglected nonlinear terms involving w c . Higher-order shear deformation theory has the 
advantage of fewer unknowns than the three-dimensional theory and yet it gives comparable results to 
those given by three-dimensional flexibility theory. For the ±45° composite plates, results show more 
pronounced nonlinear behavior in the postbuckling range as the plate width-to-thickness ratio b/h 
decreases. This more pronounced nonlinear response is a direct result of the increase in shear flexibility 
of the more compliant ±45° composite plates. The difference in the order of the approximation of the 

four theories is most evident for the ±45° composite plate results, particularly for the average axial stress 
resultant N Xav as a function of the applied displacement U. and for the maximum out-of-plane 

displacement w as a function of the applied displacement U. 

CONCLUDING REMARKS 

This paper presents buckling and postbuckling results for aluminum plates and ±45° composite plates 
subjected to longitudinal compressive end-shortening displacements. The side edges of the plates are 
simply supported and free to slide along the edges to allow constant longitudinal strain. The effects of 
varying plate width and thickness on the buckling stress coefficient is described. The buckling results for 
aluminum plates with width-to-thickness ratios less than ten, indicate that including the effects of 
transverse shear deformation is important when determining the compressive buckling stress and these 
effects should be included. Postbuckling results for plates with transverse shearing flexibility are 
compared to results from classical theory for various width-to-thickness ratios. Characteristic curves 
indicating the average longitudinal direct stress resultant as a function of the applied displacements are 
calculated based on four different theories: classical von Karman theory, a first-order shear deformation 
theory, a higher-order shear deformation theory that satisfies the bounding conditions at the upper and 
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lower surfaces of the plate, and a three-dimensional flexibility theory that can predict the transverse and 
inplane stresses directly. 

Present results indicate that the three-dimensional flexibility theory gives the lowest buckling loads for 
the four theories considered, and produces acceptable results except when the effect of the missing 
nonlinear terms involving the coefficient w c of the trigonometric term in the expansion of the transverse 
displacement w and the influence of the bounding conditions becomes dominant. The three- 
dimensional flexibility theory has the potential for permitting the development of a rigorous approach for 
obtaining direct through-the-thickness stress components without the current limitations caused by 
using additional trigonometric terms in z in the expansion of the transverse displacement w and by 
retaining the currently neglected nonlinear terms involving w°. The higher-order shear deformation 
theory has fewer unknowns than the three-dimensional flexibility theory but cannot predict transverse or 
inplane stresses. The figures presented show that, for postbuckling of aluminum plates, small 
differences occur in the average longitudinal direct stress resultant, in the maximum values of the other 
stress resultants, and in the maximum transverse displacements calculated when the effects of 
transverse shear flexibility are included in the various theories. For the ±45° composite plates, results 
show more pronounced nonlinear behavior in the postbuckling range as the plate width-to-thickness 
ratio b/h decreases. This more pronounced nonlinear response is a direct result of the increase in shear 
flexibility of the more compliant ±45° composite plates. The difference in the order of the approximation 

of the four theories is most evident for the ±45° composite plates results, particularly for the results for 
the average axial stress resultant N Xav as a function of the applied displacement U, and for the maximum 

out-of-plane displacement w as a function of the applied displacement U. 
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APPENDIX 

Governing differential equations are derived in more detail in this appendix for the two higher-order theories 
considered in this paper. 

Hioher-order shear deformation theory 


The displacements used in this theory are given by equations (3) as 


u(x,y,z) = u (x,y) + 


u 3 (x.y)-y-(y-|p- +w°(x,y)j (l) 


v(x,y,z) = v (x,y) + 


v a (x,y)-y-|^-^- +w,y(x,y)J 


(A1) 


w(x,y,z) = w°(x,y) 


Substitution of equations (3) into equation (7) gives the strain-displacement relations 

e x = u -5 + ^ w -? 2 + 

Ey = v° +iw° 2 + ~^f(v 7 + w yyh) 

Y xy = u ° + v° + w° w ° + (u a + v*) Z--± (*) 3 (u a + v,® + 2w° y h) 

rx Z f?4)(i-ef) 


(A2) 


The assumption that the displacements are sinusoidal along the length leads to 


u° = - 1)*- + uo(y) sin-2fflS- 
a * X 

u a = u a (y) cos^ 

X 

w° = w^(y) sin^L 


v° = v°(y) + vg(y) cos-2^- 

v a = v^(y) sin-®! 

X 


(A3) 



Stresses are determined from Hooke's law according to equation (8), and stress resultant forces and 
moments are determined by the following integrals through the thickness 


, x 0 + ,n x 2 


__ „ 2 Jt x 
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1- 4z 
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'-h/2 

l h / 


dz 


Nxz 0 cos^- 


y. h/2 

= I x xz 

" -h/2 
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4z 


\\ 2 i 


dz 


N X y 2 sin- 


I 7t X _ C 
* | 


l X ydZ 


sin 
A o 


h/2 

/ h/2 

a x zdz 


Mx 1 sin X = 
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h/2 

h/2 


---hoy— I dz 


-h/2 

h/2 


M yi sin-V- = 


M xy 0 cos 


/ n/<: 

OyZdZ 

-h/2 
y. h/2 

** = ( ---ha y (— 

XI 3 y \h, 

•'-hr 

k x _ r 

x J_i 


dz 


-h/2 
h/2 


TxyZdZ 


h/2 

h/2 


M ^CO s 


-| h, «nh 


dz 


(A4) 
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Substitution of the stresses and strains into the virtual work expression, equation (9), and performing the 
variation leads to the differential equations 


N v 


u 2' = ^ v 2-M P 2®- + 7^ 


X 2 A « 66 

v|' =lp 2 2 - ^(u22il + UP 2 (a-f) + ^2. 
4 ^22' ^ ^ ^ / A 22 


”^ 66 +B66 ^ P 2 ^ D 66 + 2D 66 + B 6 e ) 

± = 

h 


B 22( M y 0 _M yi ) _D 22 M y-| + ( D 12 B 22“ D 12 D 22)y^- + 


( D 1 2 B 22 " B 1 2 ° 22 ) | u 1 ^ + W 1 

A 1 2 
a 22 



-2 \ 

ID22 ^22 " ^22/ 




a 22 


N »/2 "X Nx2 

N»2-X N ^ 


M 


xy,— f-(Mx 0 +M x,)* A 55' 


a 

U 1 °7t 

TT^r 1 


M yo-f Mx »i* A «hf +P2 ) 


N y 0 ' 0 


°y’’(x) 2 ( Mx o + ^x 2 K-Mx,(f) - 


D: : 

where [>.: = - — i- in. B: : = - 
IJ 5b IJ 


D:: 


21 b 


in 2 


(AS) 


and the superscript (’) denotes differentiation with respect to y. 
The stiffnesses of the plate are given by 



(AS) 
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where the Cy are the stiffnesses in equation (8). Using the definition P 2 - w 1° ’ 9 ives * he following two 
differential equations, which complete the set of equations (14 equations with 14 unknowns) without 
squares of derivatives of the unknowns as required in the solution procedure (reference 11). 


f>2 — 


a 

v l‘ 


a 

D22(My 0 -M yi )-D 2 2My 1 + (D22Di2 -D 22 D 12)-j^-^- 


-D, 2 ,) 


(°22 D 12 ” D 22 B 1 2) "jff +W l(f') 


>22 B 22 


(A7) 


The boundary conditions at y = 0and y = bare 

u 2 = u 1 = 0 - Ny o = Ny 2 = 0, W? = 0, M yo = M yi =0 (A8) 


Three-dimen sional flexibility theory 


The displacements used tn this theory are given by equations (4) as 

u(x,y,z) = u°(x,y) + u a (x,y) L+ u s (x,y) sin 3^ 

h n 

v(x,y,z) = v°(x,y) + v a (x,y) i + v s (x,y) sin XL 

h n 

w(x,y,z) = w°(x,y) + w c (x,y) cos-^- 


Substitution of equations (4) into equations (7) and neglecting the nonlinear terms involving w° gives the 
strain-displacement relations 



(A10) 
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The assumption that the displacements are sinusoidal along the length leads to 

u 0 = -U*- + U 2 (y) sin-2nx. v° = v°(y) + v^y) cos-&K. 

3 X A, 

u a = u a (y) cos^L V a = v^(y) sin^ 

u s = u a (y) cos-B& v s = v^(y)sin-^. (All) 

X A, 

w° = w^(y) sin -EL 
X 

w c = w^(y) sin^BL 
X 


Stresses are determined from Hooke's law according to equation (8), and stress resultant forces and 
moments are determined by the following integrals through the thickness 



(A12) 
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h/2 


L Zl = 


-y*i 


Lxzi = 


-xyi 
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af sin 2 S2_dz 
z h 


xS z cos 2 ®^ dz 
yz h 


xS z cos 2a ^ dz 
xz h 


(tx a y(^)*4 slnl ^) sin Y dZ 


where the form of the stresses are 


ct y = o® + a?- 2 - +o5sin/t-2- 
x x x h h 

<V- * °yf +a y sin! ^ 


a, = cSsin®^- 
2 z h 


x yz * x yz +TyzCOS^ 

T xz = T xz +x xz 005 

T xy = T xy + ^y f + s,nB jf 


Substitution of the stresses and strains into the virtual work expression, equation (9), and performing the 
variation leads to the differential equations 

u o- 8,0 


2 .^12. 


A M 66 


h X 
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(A13) 
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N 


xy 2 


h 

= 2jln 


M4 


H> 



x 2 


N y 2 ' = _ ^ LN xy 2 
M X yi'= Nxz-| -Mx^ 

Myi' = N yzi +M xyi S. 


H 


44 


H 44 A 44 



The stiffnesses of the plate are given by 


A r}_, 

°r{' 


Ciidz 


r h/2 

^’La 


Ciicos^dz 
'J h 




h/2 

h/2 


CijZ 2 dz 


W- 


h/2 

h/2 


CiiZsinSZ-dz 
h 


h/2 

h/2 

h/2 


C::Sin S2.dz 
11 h 


(A14) 


r n 

^ * i— h/2 1 


CiiCOS 2 ®2-dz 


where the Cy are the stiffnesses in equation (8). The definition Pi = wi° ' is not used in this theory due to 
differences in the formulation of the theories. Instead, the following definition of Pi is used 



which results in the completely defined set of equations (16 equations with 16 unknowns) without squares of 
derivatives of the unknowns as required in the solution procedure. 

The boundary conditions used at y = 0 and y = bare 


u 2 = U 1 = U 1 = 0. N y =N y =0, = Wi = 0, M^-L^-0 


(A16) 
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Fig. 2 Critical axial stress resultant, N xcr , versus plate length-to-width ratio, a/b, for different 
width-to-thickness ratios, b/h, and different theories for aluminum plates of finite length. 
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Fig. 3a Average axial stress resultant, N xav , versus applied end-shortening, U, for different width-to-thickness ratios, 
b/h, and different theories for aluminum plates of infinite length where U/a is finite. 
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Fig. 3b Average axial stress resultant, N xav , versus applied end-shortening, U, for different width-to-thickness ratios, 
b/h, and different theories for ±45° composite plates of infinite length where U/a is finite. 
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Fig. 4a Maximum transverse stress resultant, Ny max . versus applied end-shortening, U, lor ditterent width-to-thickness 
ratios, b/h, and different theories lor aluminum plates of infinite length where U/a is finite. 
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Fig. 4b Maximum transverse stress resultant. Ny max , versus applied end-shortening, U, for different width-to-thickness 
ratios, b/h, and different theories for ±45° composite plates of infinite length where U/a is finite. 
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Fig. 5a Maximum shear stress resultant, N xymax , versus applied end-shortening, U, for different width-to-thickness ratios, b/h, 
and different theories for aluminum plates of infinite length where U/a is finite. 
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Fig. 5b Maximum shear stress resultant, N x y max , versus applied end-shorlening, U, for different width-to-thickness ratios, b/h, 
and different theories for ±45° composite plates of infinite length where U/a is finite. 
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Fig. 6a Out-of-plane displacement, w, versus applied end-shortening, U, for different width-to-thickness ratios, b/h, 
and different theories for aluminum plates of infinite length where U/a is finite. 
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Fig. 6b Out-of-plane displacement, w. versus applied end-shortening, U, for different width-to-thickness ratios, b/h, 
and different theories for ±45° composite plates of infinite length where U/a is finite. 
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